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Abstract 

Nonequilibrium thermodynamics has shown its applicability in a 
wide variety of different situations pertaining to fields such as physics, 
chemistry, biology, and engineering. As successful as it is, however, 
its current formulation considers only systems close to equilibrium — 
those satisfying the so-called local equilibrium hypothesis. Here we 
show that diffusion processes that occur far away from equilibrium 
can be viewed as at local equilibrium in a space that includes all the 
relevant variables in addition to the spatial coordinate. In this way, 
nonequilibrium thermodynamics can be used and the difficulties and 
ambiguities associated with the lack of a thermodynamic description 
disappear. We analyze explicitly the inertial effects in diffusion and 
outline how the main ideas can be applied to other situations. 
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Concepts of everyday use like energy, heat, and temperature have ac- 
quired a precise meaning after the development of thermodynamics. Ther- 
modynamics provides us with the basis for understanding how heat and work 
are related and with the rules that the macroscopic properties of systems at 
equilibrium follow Outside equilibrium, most of those rules do not apply 
and the aforementioned quantities cannot unambiguously be defined. There 
is, however, a natural extension of thermodynamics to systems away from 
but close to equilibrium. It is based on the local equilibrium hypothesis, 
which assumes that a system can be viewed as formed of subsystems where 
the rules of equilibrium thermodynamics apply. Due to the usual dispar- 
ity between macroscopic and microscopic scales, most systems fall into this 
category. This is the case of, for instance, the heat transfer from a flame, 
the flow through a pipe, or the electrical conduction in a wire. Nonequi- 
librium thermodynamics then extracts the general features, providing laws 
such as Fourier's, Fick's, and Ohm's, which do not depend on the detailed 
microscopic nature of the system [|]. 

In contrast, there are other situations where the local equilibrium hypoth- 
esis does not hold. Many examples are present in the relaxation of glasses 
and polymers [0,f|||, in the flow of granular media ||, and in the dynam- 
ics of colloids fl7|. The main characteristic of such systems is the similarity 
between microscopic and macroscopic scales, which usually involve internal 
variables with "slow" relaxation times. The so-called inertial effects in dif- 
fusion processes are perhaps the simplest and most illustrative example. In 
this case, the relaxation of the velocity distribution and changes in density 
occur at the same time scale. Therefore, local equilibrium is never reached. 
Here we show how nonequilibrium thermodynamics, as already established 
in the sixties H,@], can be applied to this situation. 

Nonequilibrium thermodynamics assumes that the definition of en- 
tropy S can be extended to systems close to equilibrium. Therefore, entropy 
changes are given by the Gibbs equation: 

TdS = dE + pdV - fidN , (1) 

where the thermodynamic extensive variables are the internal energy E, the 
volume V, and the number of particles N of the system. The intensive 
variables (temperature T, pressure p, and chemical potential (i) are functions 
of the extensive variables. Local equilibrium means that the Gibbs equation 
holds for a small region of the space and for changes in the variables that are 
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actually not infinitely slow. Therefore, the internal state of the system has 
to relax to equilibrium faster than variables change. In this way, all variables 
retain their usual meanings and the functional dependence between intensive 
and extensive variables is the same as in equilibrium. 

Following this approach, nonequilibrium thermodynamics has been ap- 
plied to study diffusion processes. The simplest case takes place in one di- 
mension at constant temperature, internal energy, and volume. In this case, 
from Eq. (|l|) we obtain a Gibbs equation that depends only on the density 
and the spatial coordinate x: 

Tds(x) = —fi(n,x)dn(x) . (2) 

Here s is the entropy per unit volume and n is the density. The chemical po- 
tential has the same form as in equilibrium. For instance, for an ideal system 
— formed of non-interacting particles — it is proportional to the logarithm of 
the density plus terms that do not depend on the density 0. Notice that 
these terms can include thermodynamic variables such as temperature and 
internal energy, and also the spatial coordinate. In the case of non-interacting 
Brownian particles, its explicit expression is 

li= ^L\ nn + C (x) , (3) 
m 

where m is the mass of the particles, ks the Boltzmann constant, and C(x) 
a function that takes into account possible spatial inhomogeneities. The 
dynamics of n is restricted by the mass conservation law and therefore follows 

dn dJ 

~dt = ~dx~ 1 { > 

with J being the flux of mass. An additional assumption of nonequilibrium 
thermodynamics is that this flux is given by 

where L is the phenomenological coefficient. From this, we obtain the usual 
diffusion equation 

dn _ d ( dn\ 
dt dx \ dx J 

with the diffusion coefficient D = L{dji/dn). 
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When inertial effects are present, changes in density occur at a time 
scale comparable with the time the velocities of the particles need to relax 
to equilibrium. The Gibbs equation as stated in Eq. (JJ) is no longer valid 
because local equilibrium is never reached. The entropy production depends 
also on the particular form of the velocity distribution. Both the spatial 
coordinate, x, and velocity coordinate, v, are needed to completely specify the 
state of the system. Therefore, we consider that local quantities are function 
of both coordinates. If the system is coupled to other degrees of freedom 
that relax faster than the velocity and density, a thermodynamic description 
is still possible. For instance, this is the case of Brownian particles, where 
the host fluid provides these thermodynamic degrees of freedom. Thus, we 
consider that diffusion takes place in a two-dimensional space (x, v) instead of 
in the original one-dimensional space (x). In this case, the chemical potential 
for an ideal system (e.g., non-interacting Brownian particles) is given by 

fi = lnn(x, v) + C(x, v ) , (7) 

m 

where C(x, v) is a function that does not depend on the density ||]. The 
form of this function can be obtained by realizing that at equilibrium the 
chemical potential is equal to an arbitrary constant. If we set this constant 
equal to zero, we obtain 

ksT 1 2 

u= In n -\ — v . (8) 

m 2 

Therefore, the Gibbs equation is now 

Cfcj T 1 \ 
In n(x, v) + —v 2 J dn(x, v) , (9) 

The idea of applying the rules of thermodynamics in an internal space was 
already proposed by Prigogine and Mazur and has been used in several 
situations |2|, [Ttf . In all of them, however, there was no thermodynamic 
coupling of these internal degrees of freedom with the spatial coordinate. 
This is precisely the situation we are considering here. 
In the (x, v )-space the mass conservation law is 

dn d.J x dJ v , , 

— = . (10) 

dt dx dv 
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Following the standard thermodynamic approach, the flux of mass is given 

by 

■^x — —L xx — L xv —- , (11) 
ox ov 

j -t ?t_r ^ 

where L^, with i,j = {x, v}, are the phenomenological coefficients. There 
are some restrictions on the values can take. Since the system is at local 
equilibrium in the (x, w)-space, Onsager relations imply that L xv = —L vx . 
In addition, the flux of mass in real space, J x (x) = Pf^ v n(x, v)dv, has to 
be recovered from the flux in the (x, w)-space by contracting the velocity 
coordinate: J x (x) = /f^ J x (x,v)dv. Therefore, 



/oo 
-oo 



°° a [°° ( T k B Tldn k B Tldn \ 

vndv = - L xx — + L XV — + L xv v]dv. (13) 

J -00 \ m nox m n ov I 



Since n(x, v) can take any arbitrary form, the last equality holds if and only 
if L xx = and L xv = —n. Thus, the only undetermined coefficient is L vv , 
which can depend explicitly on n, x, and v. 

Previous equations can be rewritten in a more familiar form by identifying 
the phenomenological coefficients with macroscopic quantities. In this way, 
with L vv = 11/ t, the fluxes read 

+ f!;)». (14) 

(D 8 v D d\ 

y t ox r t 2 ov J 

where D = ^-r and r are the diffusion coefficient and the velocity relaxation 
time, respectively. The equation for the density is given by 

dn d d ( v D d \ , 

dt dx dv \t t 2 dv J 

This kinetic equation is equivalent to the Fokker-Planck equation for a Brow- 
nian particle with inertia since, in an ideal system, the density is proportional 
to the probability density; i.e. n(x,v) = mNP(x,v), where P(x,v) is the 
probability density for a particle to be at x with velocity v, and N is the 
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number of particles of the system. The resulting Fokker-Planck equation 
could have also been derived by following standard techniques of stochas- 



tic processes |TT]| or kinetic theory ||12|| . which are among the microscopic 
statistical theories for studying nonequilibrium phenomena. 

The approach we have followed, however, explicitly illustrates how ther- 
modynamic concepts can be transferred from equilibrium, through local equi- 
librium, to far from equilibrium situations. The condition of equilibrium is 
characterized by the absence of dissipative fluxes (J x = and J v = 0). There- 
fore, from Eq. flliD we obtain that the velocity distribution is Gaussian with 
variance proportional to the temperature. If deviations from equilibrium are 
small (J x 7^ and J v — 0), the local equilibrium hypothesis holds. This is 
the domain of validity of Fick's law, 

J. - , (IT) 

which is obtained directly from the equations for the fluxes. In this case, 
the distribution of velocities is still Gaussian but now centered at a non- 
zero average velocity and the variance of the distribution is related to the 
temperature in the same way as in equilibrium. Beyond local equilibrium 
( J x 7^ and J v ^ 0), the velocity distribution can take any arbitrary form, 
from which there is no clear way to assign a temperature. There is, however, 
a well defined temperature T: that of local equilibrium in the (x, v )-space. 

In Fig. [I] we illustrate the concepts discussed previously. We show the 
velocity profiles obtained from Eq. (^) for two representative situations. 
For fast relaxation of the velocity coordinate, the velocity distribution is 
Gaussian and centered slightly away from zero, in accordance with local 
equilibrium concepts. For slow relaxation, however, the velocity distribution 
loses its symmetry (and its Gaussian form). In this case, the temperature 
does not give directly the form of the distribution and one has to resort to 
local equilibrium in the (x, v )-space to describe the system. 

It is important to emphasize that the temperature T is the one that 
enters the total entropy changes and therefore the one related to the second 
principle of thermodynamics. Other definitions of temperature are possible 
though. To illustrate this point, let us compute the entropy production a. 
This quantity is obtained from local changes in entropy, which are given not 
only by the production but also by the flow: 

rr^s dn ( dJ Sx dJ Sv \ . . 
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where {Jsx, Jsv) is the entropy flux. In our case, the expression for the 
entropy production is 



n(x,v) ( k B Td\nn(x,v)\ 

aM = -Fr{ v + — dv ) ■ (19) 

Now, given a Gaussian velocity distribution n(x,v) = n (x)e~ mv2 ^ 2kBT( - x \ we 
can easily understand the meaning of the temperature T(x) defined through 
the variance of the distribution: it is the temperature at which the system 
would be at equilibrium (a = 0). The definition of an effective temperature 
as that giving zero entropy production can be extended to arbitrary velocity 
distributions. From Eq. (|l9|) , we obtain 

i 1 ( d In n(x,v)\ , . 

k B • (20) 



f(x,v) vm \ dv 

The temperature defined in this way is formally analogous to the equilibrium 
temperature since the right hand side term of the preceding equation can be 
rewritten as the derivative of an entropy with respect to an energy: 



dsJx, v) 



T(x,v) de(v) 



(21) 



where s c (x,v) = —^\nn(x,v) and e(v) = \v 2 . The term s c (x,v) and e(v) 
can be viewed as the configurational entropy and the kinetic energy per unit 
of mass, respectively. In general, other definitions of effective temperature are 
possible. For instance, by considering e(v — v(x)) instead of e(v) in Eq. (pT|), 
the resulting temperature would be that of local equilibrium. In this case, 
however, this temperature does not give zero entropy production but just 
that of the macroscopic motion. This temperature is then the one at which, 
once the macroscopic motion is disregarded, the internal configuration of the 
system would be at equilibrium. 

In general, since T(x, v) is not only a function of x but also of v , given 
a point in space, there is no temperature at which the system would be at 
equilibrium; i.e. T(x, v) ^ T(x). If an effective temperature at a point x were 
defined, it would depend on the way the additional coordinate is eliminated. 
Thus, ambiguities in far from equilibrium quantities arise when considering 
a lower dimensional space than the one in which the process is actually 
occurring. This is to some extent similar to what happens with effective 
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temperatures defined through fluctuation-dissipation theorems. In such a 
case, the effective temperature can depend on the scale of observation |i~3|| . 
It is interesting to point out that all these effective temperatures, despite 
their possible analogies with the equilibrium temperature, do not have to 
follow the usual thermodynamic rules since the system is not actually at 
equilibrium at the temperature T. 

The idea of increasing the dimensionality of the space were diffusion takes 
place, so to include as many dimensions as non-equilibrated degrees of free- 
dom the system has, can also be applied to other situations. In a general 
case, the additional degrees of freedom do not necessarily correspond to the 
velocity. For instance, let us consider a degree of freedom Q(x) that at local 
equilibrium enters the Gibbs equation in the following way: 

Tds(x) = -n dn{x) - B dG(x) , (22) 

where B = B(n, 6, T) = T(ds/dQ) nt T- In this case, one usually assumes that 
given T , n(x), and Q(x), the function B is completely determined through 
the equilibrium properties of the system. Far away from equilibrium, we 
would have to consider explicitly an additional coordinate 9, which is related 
to the degree of freedom by Q(x) = J 9n(x,9)d9. The corresponding Gibbs 
equation 

Tds(x,9) = —fidn(x,8) (23) 

would have to take into account the dependence on the coordinate 8 through 
the chemical potential //. Once the Gibbs equation has been obtained, the 
way to proceed would be analogous to the one we followed for the inertial 
effects. For instance, some systems with both translational and orientational 
degrees of freedom can be described by the chemical potential 

m f{6) 

where 6 is now an angular coordinate, U cos 6 is the orientational energy, 
and f{6) is a function accounting for the degeneracy of the orientational 
states (for rotation in three and two dimensions, f{6) = sinO and f{6) = 1, 
respectively) |2|]. This type of systems include, among others, liquid crystals 
and suspensions of rod-like particles ||, field-responsive suspensions PH| , and 
polarized systems ||. At local equilibrium some instances of B and 9 are 
then electric field and polarization; and magnetic field and magnetization. 



Beyond local equilibrium, by writing the (x, 6) counterpart of Eqs. (|10|) . 



8 



(|TT|), and flT2]), one can obtain a kinetic equation that describes the behavior 
of the system. This equation includes as particular cases the Fokker-Plank 
equations obtained for those systems by means of microscopic theories ||o|,|To 



Along this paper, we have been assuming ideality and locality. The condi- 
tion of ideality is that the system consists of non-interacting particles. In this 
case, the chemical potential is proportional to the logarithm of the density 
plus terms that do not depend on this quantity. Non-ideality can be directly 
taken into account by considering the right dependence of the thermodynamic 
quantities on the density and, in general, will give rise to nonlinear partial 
differential equations. A more difficult aspect to deal with is the presence of 
non-local effects. In such the interactions between the different parts 

of the system will need of integro-differential equations to be incorporated in 
the description. 

The main result of our analysis shows that in far from equilibrium dif- 
fusion processes, local equilibrium can be recovered when all the relevant 
degrees of freedom are considered at the same level as the spatial coordi- 
nate. In the resulting extended space, thermodynamic quantities, such as 
temperature and the chemical potential, admit a well defined interpretation. 
The scheme we have developed may then provide the basis for a consistent 
formulation of thermodynamics far from equilibrium. 
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Figure 1: Velocity profiles obtained from Eq. (|i6|) when a density gradient 
is applied. The solution has been obtained through a standard numerical 



algorithm following a first order upwind discretization scheme ||16|| . The 
system is in a rectangular domain in the (x, v )-space, from x = to x = 
1, and from v — — 10 to f= 10. The lower and upper dashed curves in 
the figure represent the boundary conditions applied at x = and x = 
1, respectively: n(l, v) = 10n(0,t>) = (27r)~ - 5 exp(— v 2 /2). Filled circles 
correspond to velocity profiles at x = 0.5 for fast relaxation of the velocity 
coordinate (r = 0.1), whereas empty circles correspond to slow relaxation 
(r = 10). In both cases D/r = ksT /m = 1. All values are given in arbitrary 
units. 
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